Destruction of interference by many-body interactions in cold atomic Bose gases 
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We study the effects of many-body interactions on the interference in a Mach-Zehnder interfer- 
ometer for matter waves of ultracold Bose atoms. After switching off an axial trapping potential, 
the thermal initial wavepacket expands, and subsequently interference fringes may be observed in 
a circular ID trap. These are computed for axial harmonic or 5-function traps, and for interaction 
strengths from the Thomas-Fermi regime to the Tonks- Girardeau limit. It is shown that many- 
body correlations in a realistic setup destroy interference to a large degree. Analytical expressions 
allowing to infer the observability of phase coherence and interference are provided. 
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I. INTRODUCTION 

Recent experimental advances in the control of ultra- 
cold atomic gases have opened new exciting possibilities 
to systematically study the effects of many-body interac- 
tions in low-dimensional strongly correlated mesoscopic 
systems Q- 0- ■ Two main experimental routes are 
presently being pursued, namely magnetic or magneto- 
optical trap technology @, |j, [a, Q ,and the formation of 
atomic waveguides on microchips 0, H S LLi3 • A partic- 
ularly interesting limit arises in the true ID case, which 
may, for instance, be realized by confining N atoms in a 
highly anisotropic cigar-shaped trap with large transverse 
trapping frequency u>± [TIL IT^| . In this paper we study 
interacting ID Bose gases, but similar ideas and conclu- 
sions apply to cold fermionic gases as well. In the ID case, 
in an infinitely long system, quantum fluctuations pre- 
vent the existence of a Bose-Einstein condensate (BEC) 
even at T = 0. For a finite system and very weak in- 
teractions, however, BEC could happen at extremely low 
temperatures not reachable at present |Iall4llalldll7|| . 
The physics of a trapped ID atomic Bose gas is then char- 
acterized by a crossover from the high-density weakly 
correlated Thomas-Fermi (TF) regime to a strongly cor- 
related Tonks-Girardeau (TG) regime 0,0] at low den- 
sity, see Refs. 0,123, El EJ. The relevant dimensionless 
interaction parameter governing this crossover in a uni- 
form system is 

7 = mg/h 2 p , (1) 

where g is an effective ID interaction constant [2j|, see 
Eq. © below, and po is the density. While for the 
uniform system, a Bethe ansatz solution due to Licb 
and Liniger [24| is available, for the case of a (possibly 
time-dependent) axial trapping potential V{x,t), no ex- 
act statements are known, and one generally resorts to 
approximate or numerical methods 25, 26U27|. We men- 
tion in passing that this system also has very interesting 
single-particle excitations with exotic properties such as 
fractional (anyon) statistics. While such systems could 
also be realized in more conventional condensed matter 
systems, the high quality and tunability of atomic Bose 



gases renders them very attractive for fundamental stud- 
ies. 

In this paper, we address the influence of interactions 
on the phase coherence and the interference properties 
of a ID Bose gas. To be specific, we consider the cir- 
cular geometry indicated in Fig. ^ At time t = 0, the 
axial trapping potential V(x,t) localized around x « 
is switched off, and after some characteristic expansion 
time, the expanding right- and left-moving wave packets 
will then meet at the other side of the ring, x « ±L/2. 
If phase coherence of the wave packets is not lost during 
the propagation, a Mach-Zehnder-type interference sig- 
nal should be observable at this side. While a modified 
Gross-Pitaevskii (GP) approach, obtained from combin- 
ing the Lieb-Liniger solution with a local density approx- 
imation (LDA) valid for large enough systems [2(j , is able 
to yield accurate ground-state density profiles po(%, t) of 
a ID Bose gas after switching off the axial trap even 
in the Tonks limit [2{|, it typically overestimates in- 
terference signals by orders of magnitude j2(|. For a 
correct description of interference, it is thus necessary 
to take into account quantum-mechanical phase fluctua- 
tions, which have so far only been addressed in the static 
case pgl l29j or in the Tonks limit under rather special 
initial conditions [sol Isij . Below we show that under 
realistic modelling of the initial preparation and the sub- 
sequent dynamics, phase coherence and interference are 
strongly reduced by many-body interactions in ID Bose 
gases. Moreover, we provide quantitative expressions to 
assess this effect in experiments. Proper insight into this 
question is of practical importance for the development 
of future atom interferometers and for the understand- 
ing of current experiments [lfil l32l . The remainder 
of the paper is organized as follows. In Sec. [n] we study 
magneto-optical or magnetic traps with 

V(x,t) = ~mu 2 x x 2 G(-t), (2) 

with angular frequencies lo x uj±; here O is the Heavi- 
side function. In Sec.|nl the validity of LDA is presumed. 
In Sec. IIIII we briefly discuss a related setup employ- 
ing a (5-trap which may be more realistic for microchip- 
trapped gases yet allows for an exact solution of the time- 
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FIG. 1: Schematic geometry of a Mach-Zehnder interfer- 
ometer realized in a circular trap with circumference L. By 
applying an axial trap potential at t < around x ~ 0, a 
thermal non-equilibrium initial state is prepared, whose time 
evolution and interference at x ~ ±L/2 are the subject of this 
paper. 



dependent problem. Finally, we conclude in Sec. IIVI 



II. INTERFERENCE SIGNAL: HARMONIC 
INITIAL TRAP 

Let us consider a circular trap at ultralow temperature, 
ksT <C fiwi, such that motion in the transverse direction 
is frozen out, see Fig. ^ The effective ID interaction 
strength is [23| 



2h 2 



a 1D = --±[l-C(a s /d ± )], (3) 
tociid za s 

where a s > is the 3D s-wave scattering length de- 
scribing the repulsive atom-atom interactions, d± = 
(2/i/mt^) 1 / 2 , and C w 1.4603. The interaction strength 
can be tuned by Feshbach resonances [34[ or by using 
optical lattices ■ 



A. Fluctuation modes 

The resulting Hamiltonian of the ID Bose gas at 
fixed particle number AT, including an arbitrary time- 
dependent axial trap potential V(x,t), is 



H(t) 



L/2 



dx $ (x, t) 



-L/2 

— I dx 



2m 



d 2 x + V{x,t) 



4>(x,t) 

(4) 



where the colons denote normal-ordering. For V(x, t) = 
0, this model was solved exactly by Lieb and Liniger [2J|. 
To make progress in the presence of the external poten- 
tial, the field operator ip is expressed in terms of the 
density p(x) and the phase </>(x), 



%j)(x) = yj p(x)t 



(5) 



We employ canonically conjugate density [phase] fluc- 
tuation operators Ii{x,t) — p(x,t) — p (x,t) [Q(x,t) — 
4>(x, t) — (f>o(x, t)] describing quantum fluctuations around 
the solution ipo(x,t) = \/ po(x, t) exx>\i(j)c\(x. t)] of the 
time-dependent modified GP equation 20, 25], 



ihd t ^[){x,t) = 



-^d 2 x + V(x,t) + F(p (x,t)) 



(6) 

where F(po) is the (exactly known) chemical potential 
for the Lieb-Liniger problem |2f| . with limiting values 



F(p ) = 



5Po, 
TT 2 h 2 pl/2m, 



7 < 1 
7 > 1. 



(7) 



For arbitrary time-dependent potential V(x,t), no closed 
solution of Eq. © for arbitrary 7 is available. From nu- 
merical work, however, the expansion of a ID Bose gas 
after switching off an harmonic axial trap potential has 
recently been shown to violate self-similarity |25| in the 
intermediate regime between the TF and the TG limit, 
7 w 1, which can be traced back to the variable nonlin- 
earity of F(po) for different 7. 

Substituting Eq. J5J into Eq. J3J and expanding up 
to quadratic order in $ and II, we obtain the effective 
low-energy Hamiltonian describing quantum fluctuations 
around the solution to Eq. ijfjj . Neglecting terms contain- 
ing n<9^n or (p " 1 9^po)n 2 , we obtain the time-dependent 
Hamiltonian 



H(t) 



dx 



h 2 Po 



2m 



(8) 



m 2 opo 



While terms involving n9 2 n are irrelevant in the long- 
wavelength low-energy regime, the neglect of terms ex 
(Pq d 2 po)H 2 req uires that the interaction strength g is 
not too weak [35j. In practice, this approximation is valid 
to high accuracy even in the TF regime, and therefore in 
all cases of interest here. Interference fringes are correctly 
captured in such a long-wavelength theory despite their 
rapidly oscillating character, since the oscillations arise 
due to a mixing of left- and right-moving wave-packets. 
These wave-packets are both correctly described within 
the present theory. 

Equation JSJl resembles the Hamiltonian of a Luttinger 
liquid 13] , generalized to a time-dependent and spatially 
non-uniform system, where a non-standard coupling be- 
tween n and d x $ is induced by the time dependence of 
the external potential. The Hamiltonian is quadratic 
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in the fluctuation operators and can therefore be diago- 
nalized by a Bogoliubov-de Gennes transformation. How- 
ever, explicit diagonalization is rather difficult for arbi- 
trary potential V(x,t), and here we shall limit ourselves 
to the form given in Eq. ©. Although the full Hamilto- 
nian is time- independent for t > 0, the Hamiltonian 
governing quantum fluctuations becomes explicitly 
time-dependent due to the time dependence of po(x,t) 
and 4>o{x, t). Furthermore, to make progress, we restrict 
ourselves to the self-similar limits 7 1 and 7 ^> 1. 
Then the densit y P r\(x,t > 0) is self-similar with some 
scale factor b(t) |25| . 



p (x,t) 



n(x/b(t)) 



(9) 



where b(0) = 1 and n(x) — po(x,0) is the initial density 
profile. Solving Eq. 10 with Eq. © results in the phase 



<MM) 



mx 2 b(t) 
~~2h~b(t)' 



(10) 



Under self-similarity, it is then possible to diagonalize 
the quadratic Hamiltonian (JHJ) by solving the equation of 
motion generated by it. Using the long derivative 

D = d t +x(b/b)d x , (11) 

some algebra shows that the equation of motion for $ is 



D 



al\m(x,t) = -^^-d x [p d x ^], (12) 
b I m po 



where a = 1 in the TF regime and a = 2 in the TG 
regime. To solve Eq. (|12|) . we use the ansatz 

=X) cr i(*)/i(«/6(*))4-(*) +H.c, (13) 



where B] , Bj are standard Bose operators with time 
dependence governed by time-dependent eigenfrcqucncy 



Bj(t) = cxp 



-1 I dt'ttjit') 
'0 



B 3 (0). 



(14) 



Furthermore, fj are suitably normalized eigenfunctions, 
and the prefactor Cj is needed to bring H(t) into the 
canonical form 



H(t)=J2^ 1 (t)B}(t)B,(t), 



(15) 



where the summation extends over all eigenmodes. One 
finds a similar expression to Eq. (|13|) for the conjugate 
field n by using the commutation relation with $. Given 
the scale-factor b(t) and the initial density n(x) enter- 
ing Eq. 10, solving the eigenproblem (|12fl thus leads to 
a complete description of the time-dependent quantum 
phase fluctuations. This program is carried out in this 



section separately for the TF limit, see Sec. Ill Bl and the 
TG regime, see Sec-UHl 

To investigate coherence properties during the self- 
similar expansion process, we calculate the (equal-time) 
density matrix 



W{x,x',t) = (ft(x,t)ip(x',t)). 



(16) 



Following Refs. 0, |3(|, here it is justified to neglect 
density fluctations against phase fluctuations, resulting 
in 

W(x,x',t) ~ W {x,x' 1 t)e- F{x ^ x '' t \ (17) 
W (x,x',t) = ^po(x 1 t)p (x',t)e-^' t ^ x '^\ 
F(x,x',t) = ([<f>(x,t) - §{x' ,t)] 2 )/2. 

The task is then to compute the fluctuation correlator F. 
Substituting Eq. (|13|) into F(x,x',t), it follows that 

oc 

F(x,x',t) = Y J ^\CAt)\ 2 {l + ^n 3 )[} 3 {x,t)~f ] {x l ,t)]\ 

OO 

= \ C i (*)l 2 coth[hn 3 (t)/2k B T} 

j=0 



(18) 



x ^(x^-f.ix^t)] 2 , 

with the Bose-Einstein distribution n 3 = (B^Bj) = 
[expihQji^/kBT) -I]- 1 . 

To analyze interference fringes, it is convenient to 
switch to y = x — L/2 (y = x + L/2) for L/2 > x > 
(-L/2 < x < 0), such that the center of the overlapping 
clouds is at y = 0, and the detection signal is given by 

I(y,t) = (p(y + L/2,t)) + (p(-y + L/2,t)) 
+ 2ReW{y + L/2,-y + L/2,t). 

The last term corresponds to the (first-order) interference 
signal /int(y) of interest here. Using Eq. H17J1 . it reads 



-fint(y) = 2cos[2Qy]c 



-F(y+L/2.-y+L/2,t) 



(19) 



x y/po(y + L/2,t)po(-y + L/2,t), 



with time-dependent modulation wavevector Q — [4>a(y+ 
L/2,t) - <j> (-y + L/2, t)]/2y, which for \y\ < L takes the 
form 



Q{t) = d x <t> Q {L/2,t) = 1 ^ h l . 



(20) 



It is obvious from Eq. (|19|l that phase fluctuations enter- 
ing F(x,x',t) suppress interference. 



B. Thomas-Fermi regime 

In this subsection, we assume that the initial trapped 
gas (at t < 0) is within the TF regime, and also stays 
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within this regime throughout all relevant timescales. 
Then the initial density is given by |37| 



n(x) 



2g 



w 2 x (R 2 - x 2 )G(R 2 - x 2 ), 



(21) 



where the TF radius R follows from J n{x)dx = N, 

R = (3A^/|a 1D |) 1/3 , l x = y/h/mLj x . (22) 

Provided the expanding gas indeed stays in the TF 
regime, the expansion is well described b y th e self-similar 
profile © , where b(t) is the solution to [33, |33 



tol/b 2 , with b{0) = 1, 6(0) = 0. 



(23) 



The phase <po(x,t) is then given by Eq. (|10|l . While 
the eigenproblem for the quantum fluctuations has previ- 
ously been solved in the static TF case p^lMI^IM!!^. 
the coupling of d x $ and II in the time-dependent case, see 
Eq. |JSJ, complicates the analysis considerably. Remark- 
ably, it is nevertheless possible to solve the full dynamical 
problem analytically. 

The equation of motion l|12|) is solved under the ansatz 
JEJ with 



Cj(t) = [ZhRtt^b^/g]- 1 / 2 . 



(24) 



To see this, we note that the long derivative l|ll|) acts 
on the part belonging to the annihilation operator Bj in 
Eq. O as 



D<f>j = -m, 3 - 

J \ 1 2b 20, 



4\ 



After some algebra, using y — x/Rb(t), Eq. 112f) then 
takes the form 



d 

dy 



= o, 



(25) 



where Gj (t) is given by 



Gj = b [Qj — ^(fij/%) 



+ fljfiSlj - (b/2b) 2 + b/2b 



(26) 

Equation i|25|) has a standard solution in terms of Legen- 
dre polynomials Pj, with normalized eigenfunctions 



f j (y)= V7 + T72 Pjiv), J dyF 2 (y) = 1 , 

iff Gj (t) is time- independent and given by 
G,-[%(i)]=j(j + l)<4 



(27) 



(28) 



where j is a non-negative integer. Eigcnfrcqucncics Qj(t) 
thus follow as solutions to the differential equation in- 
duced by Eq. J5SJ). On not too long time scales, Eq. i(2*5|) 
with 1)26(1 has the solution 



ftj(*) _ / 2j(j + l) + l 
w a V 463 (*) ' 



(29) 



which holds as long as b(t)/2b(t) <C £lj(t). For the results 
shown below in Figs. [21 and [3] this condition is accurately 
fulfilled. On much longer time scales, one should instead 
(numerically) solve Eq. (|28|) . The fluctuation correlations 
(|17fl and the interference signal (|19|) then follow immedi- 
ately using these results. It is obvious that Cj(t) oc g 1 ^ 2 , 
and therefore the interaction parameter g enters the fluc- 
tuation correlator F and suppresses the interference in 
an exponential manner. The validity of the self-similar 
TF solution and the subsequent treatment of the quan- 
tum fluctuation modes is guaranteed by the validity of 
i) = 9P0i which in turn is accurate to 1% for 7 < 0.05 



To give concrete examples, consider a system of N = 
10 3 23 Na atoms on a ring of circumference L — 16R. 
where R is the TF radius. The axial trap potential is 
switched off at t = 0, and we study the resulting inter- 
ference signal. A qualitative estimate for the effective 
interaction strength at time t > can be given in terms 
of 7 (t) = 2/[\a 1D \p m (t)} [iJm, where p m {t) = p (x = 
0,t), Let us first consider a system with trap frequencies 
lu x = 0.5 kHz and co± = 50 kHz, where 7(0) w 0.0006 
indicates that the system is deeply in the TF regime and 
phase fluctuations are very small. The expanding wave 
packets begin to meet at time t = 13 ms. We calcu- 
lated the phase fluctuations from Eq. (|18fl with a UV 
cutoff hflj < fi given by the t = TF chemical poten- 
tial, fj. = mujlR 2 /2, see Eq. (JS2J|. In Fig.[2 the resulting 
interference fringes at time t = 16 ms corresponding to 
b{t) — 10 are shown. While there is a detectable inter- 
ference signal at T = 4 nK, interference is completely 
washed out by thermal fluctuations at 10 nK. As second 
example, consider the same system but with higher trap 
frequencies, lo x — 1 kHz and ui± = 100 kHz, see Fig. [3 
Here we have stronger interactions, 7(0) w 0.001, but the 
system still stays in the TF regime on the timescales of 
interest. Clearly, pronounced interference patterns can 
be observed again. 



C. Tonks-Girardeau regime 

When the system leaves the TF regime during the ex- 
pansion process, self-similarity is violated |25|. and the 
solution of Eq. 1|12[> is much more complicated. However, 
analytical progress is possible in the opposite TG limit. 
Assuming that one starts with an initial trap correspond- 
ing to the TG regime, during the subsequent expansion 
one will also stay in this regime, and the density becomes 
self-similar again [2(j, see Eq. © with b(t) = ^/l+uj^t 2 
and 



(30) 



n{x) = n° Ton ^ I - x 2 / Rl, 
where n^ onks = 2N/kR and 

R Q = v / 2Nh/muj x (31) 

denotes the Tonks radius. Furthermore, Eq. (|10fl specifies 
0o (x, t) again. 
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FIG. 2: The scaled density profile in units of p m (0) for N = 
10 3 expanding Na atoms at time t = 16 ms and temperature 
T = 4 nK. Trap frequencies are u x = 0.5 kHz and = 
50 kHz. The inset shows the corresponding density profile at 
T = 10 nK. 



0.13 




FIG. 3: Same as Fig. [5] but with oj x = 1 kHz and = 
100 kHz. The expanding clouds meet at time t = 6.5 ms, the 
snapshot is shown for t = 8 ms. The inset shows details of 
the interference fringes. 



Now Eq. i|12|) can be solved by the ansatz (|13|) with 

C 3 (t) = pfimJZofyW&W^yLto] -172 - (32) 

Proceeding in the same spirit as above, using y — 
x/Rob(t), Eq. (12) leads to 

+ (1 - y 2 )-^/^) ~ !/-£-/;(») = 0, (33) 



0.2 




FIG. 4: The scaled density profile in units of n Xonks for N = 
10 3 expanding 87 Rb atoms at time t = 1 s (see text). The 
inset shows details of the interference fringes. 

where 

Gj = b* m| - ^(hj/Slj) 2 + <">, 2L>, + 6/6^J . (34) 

With j denoting positive integers, Eq. H33fl has a solution 
provided that 

Gj = j 2 ul (35) 

which gives a differential equation determining £lj{t). 
The solution to Eq. I|33|) is then 

where P- are Jacobi polynomials, and the nor- 

malization condition is 

/ 1 i dy(l- 2 / 2 )- 1/2 /|(y) = i- 

Amazingly, Eq. (|35|l is solved by the simple form 

n j (t)/u; x = j/b 2 (t), (37) 

which holds on all timescales. It is then straightforward 
to obtain the fluctuation correlations from Eq. i|18|) . and 
the interference signal from Eq. (|19fl . 

Notably, even at zero temperature, interference is prac- 
tically completely suppressed by many-body interactions. 
In Fig. 0] we show the results of the above calcula- 
tion for N = 10 3 87 Rb atoms on a ring of circumfer- 
ence L — 16i?o, where we assume a scattering length 
a s = 42.32 nm, e.g. resulting from tuning via a Fesh- 
bach resonance. Trap frequencies are lo x — 10 Hz and 
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lu±_ = 100 kHz, which puts the system well into the 
Tonks regime, 7(0) ps 14.3. At t — 0.8 s the expand- 
ing wavepackets begin to meet, and the scaled density 
profile at time t — 1 s and zero temperature is shown in 
Fig. 0] Interference amplitudes have decreased by orders 
of magnitude when compared to the TF regime, but still 
tiny signatures are observable. At finite temperature, one 
then finds virtually no sign of interference anymore, in ac- 
cordance with expectations based on previous work [28| . 
As a result, phase coherence is drastically suppressed in 
the Tonks limit. 

Finally, in between the TF and TG limits, one must 
resort to numerical techniques. After numerical solu- 
tion of Eq. 0, numerical data for po{x,t) and 4>o(x,t) 
could be used to solve the eigenproblcm corresponding 
to Eq. (|12|) . and subsequently the interference problem. 
It is presently unclear whether this program can be car- 
ried out in practice. We shall pursue a different route 
in the next section, where analytical results are given 
for a related setup employing a d potential for the ini- 
tial axial trap. This allows to access the full crossover 
between the two limiting cases without approximations. 
These results qualitatively confirm our above statements 
obtained in the TF and TG limits. 



III. ALTERNATIVE APPROACH 

So far we have studied interference properties of atomic 
Bose gases in ID circular traps after switching off an ax- 
ial harmonic trap potential. Our treatment had to rely 
on the local density approximation (LDA) [2(j (which 
requires sufficiently small oj x ), and, moreover, explicit 
results were only obtained in the TF and TG limits. To 
check that the above picture is consistent and robust, 
in this section we consider a trap potential of the form 
V(x,t) — — VqS(x)Q(— t), assuming the circular atom 
waveguide in Fig. ^ to already contain No atoms in the 
absence of V(x, t), so that po = Nq/L is the density 
away from the trap at t < 0. Let us then imagine that 
N <C N additional atoms are injected (N oc Vq is tuned 
by the trap depth) in the distant past at x — by adi- 
abatically switching on V(x,t), e.g. from a BEC atom 
reservoir using a quantum tweezer |4l| . At time t = 0, 
the trap is switched off, and we again wish to compute 
the interference signal. This setup may be of relevance to 
experiments using injection of particles with a quantum 
tweezer, or for microchip-trapped cold atoms where con- 
finement potentials are quite steep Importantly, un- 
der the condition N <C No, this interference problem can 
be solved exactly by virtue of the bosonization method. 
Using Eq. (JSJ) and the quantum fluctuations defined in 
Sec.|n] one now arrives at a standard (uniform and time- 
independent) Luttinger liquid Hamiltonian [l3| plus the 
contribution of the trap, 

H (t) = 7^ J dxiKid^+K-^nU^-VoQi-m^t), 

(38) 



where u is the sound velocity and 1 < K < 00 is 
the dimensionless Luttinger parameter. Both parame- 
ters are determined by g and po [l3j |. In particular, in 
the TG limit, K = 1, while K > 1 in the TF limit. 
The effect of the trap potential can then be included 
exactly by the time-dependent unitary transformation 
U(x,t) -> Il(x,t) + N6(x)Q(-t), where N = KVq/ttHu 
(N is assumed to be integer) gives the number of added 
atoms. This allows us to obtain the expansion dynamics 
Po{x,t) = po + (n(a;,t)) and the interference signal in 
closed form. From the equation of motion, 



(-d'f + u 2 d 2 x )Tl{x,t) = NQ{-t)dl§{x), 



(39) 



for t < we infer the density nix) — po + N8(x), while 
for t > the expanding density profile is 

po(x, t)=p + N[S(x - ut) + S(x + ut)]/2. (40) 

Apparently, the initial ^-peak in the density is split into 
two counterpropagating parts, each moving with velocity 
u. Furthermore, the phase 4>o{x,t) is easily found in the 
form 

ttN 

The interference signal then follows from Eq. I|19(l , with 
the quantum fluctuations at x « x' w ±L/2 entering as 



Flx.x'.t) = In 

V ' ' ' 2K 



hu . , fnLk B T\ 

sinn — 

TTaoksT \ hu J 



, (41) 



where ao is a nonuniversal UV cutoff length limiting the 
applicability of the continuum model (|38|l . Note that 
due to our assumption of a constant density po in the 
absence of the axial trap, the fluctuation factor is now 
time-independent. Equation (|41|) implies power-law sup- 
pression of the interference signal oc 1/L 1 / 2K for short 
circumference L < Lt, where 



(42) 



This power law decays quite rapidly in the Tonks limit 
(K = 1), but interference is only weakly suppressed in 
the TF limit. Moreover, for large circumference L > Lt, 
we find an exponential suppression of the interference 
signal, oc exp(—L/2KLT)- The findings of the previous 
section are therefore confirmed from this simple yet exact 
calculation. There is a smooth crossover between the TF 
and TG limits, described by the above expressions. 



IV. CONCLUSIONS 

In this paper, interference properties of interacting 
atomic Bose gases in a ID circular trap have been an- 
alyzed. After switching off a trapping potential keeping 
the atoms in the narrow initial region, the expanding 
atom clouds meet at the opposite side and may produce 
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an interference signal. This signal has been computed 
for the case of a harmonic initial trap, with explicit re- 
sults in the Tonks-Girardeau regime and in the Thomas- 
Fermi limit, including the effects of thermal fluctuations. 
Our central conclusion is that many-body interactions 
strongly suppress phase coherence and interference. In 
fact, using this setup interference is unlikely to be ob- 
served in the strongly interacting Tonks gas. Even in 
the Thomas-Fermi regime, one needs rather low temper- 
atures and steep initial traps to ensure a significant inter- 
ference signal. Further support for these conclusions can 
be drawn from a simple yet exact calculation based on a 



slightly different setup where one injects a small number 
of additional atoms into the system, traps them for some 
time and then studies the expansion and the subsequent 
interference signal. 
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